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A procedure for extracting information from a heart sound signal. 

5 The invention relates to a procedure for extracting information from a 

phonocardiographic signal obtained from a transducer and subjected to signal 
processing in order to aid evaluation and diagnosis of heart conditions. The invention 
fiirthermore relates to techniques forming part of such extraction and apparatus to 
perform such feature extraction as well as coding the features to aid the ability to 

10 distinguish between related features. 

Signals obtained by means of a transducer are phonocardiographic representations of 
sounds traditionally listened to by means of a stethoscope. Training in auscultation 
takes a long time and requires an aptitude for recognising and classifying avural cues, 
15 frequently in a noisy environment. 20-30 different conditions may need to be 

differentiated, and within each, the severity evaluated. Furthermore, there may be 
combinations among these. These factors contribute to explaining why not all 
physicians perform equally well when diagnosing heart conditions, and why it may 
be time-consuming. 

20 

The so-called first (SI) and second (S2) heart sound are very important markers in 
the assessment of a heart soimd signal. These sounds are directly related to the 
fiinctioning of the heart valves, in that SI is caused by the closure of the 
atrioventricular valves and contraction of the ventricles and S2 is caused by the 
25 closure of the aortic and pulmonary valves. 

A number of patents relate to the extraction of the SI and S2 signals, such as US 
6,048,319, which concerns the measurement of the time interval between the SI and 
S2 signals in relation to the heart rate in order to determine the degree of coronary 
30 artery disease. The measurement is based on peak detection and autocorrelation and 
it may be considered a relatively slow process. 

According to the present invention the detection of S 1 and S2 is obtained by 
performing the steps of feature extraction and classification based on the energy 
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distribution over time in a feature time function. The feature extraction is performed 
by the steps of bandpass filtering, followed by instantaneous power and lowpass 
filtering. This generates a series of signal peaks or "hills", each relating to either an 
SI or an S2, and a signal classification step determines which "hill" is to be regarded 
5 as either an SI or an S2, whereby a systole is correctly identified. 



A different category of signals related to various heart conditions is generally known 
as murmurs. The known procedures of isolating and categorizing murmurs are 
generally dependent on the simultaneous recording of electrocardiographic data, such 
10 as US 5,957,866 and US 6,050,950 and this complicates the practical use of 
auscultation techniques considerably. 

The above solutions are very complex and rely on techniques that are equivalent to a 
long averaging time. According to the invention a method has been derived which is 

15 more precise and obtains a faster result. This is obtained by a sequence of steps, 
comprising an optional adaptive noise reduction, detection of SI and S2, e.g. by 
means of the feature extraction procedure mentioned above, enhancement of the 
signal by elimination of the SI and S2 contributions, performing spectral analysis 
and feature enhancement in order to obtain the energy content present in areas of a 

20 time-fi'equency representation delimited by firequency band times time interval in the 
form of energy distributions, classifying the energy distributions according to pre- 
defined criteria, and comparing the energy distributions to a catalogue of 
distributions related to known medical conditions and extracting information by 
comparing the enhanced signal to stored time functions. 

25 

The.correct.placementin.time of SI and.S2-pennitsJthe energy-xelating-to these 
sounds to be eliminated in the signal processing, and the resulting sound (including 
murmurs, etc.) is a useful starting signal for further analysis, because it increases the 
dynamic range of the remaining signal. It also permits presenting the remaining 
30 signal to the ears with a superposition of correctly placed but "neutral" SI and S2 
contributions as mere time markers, but without any signal that the listener needs to 
process in the listening process. 
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Diagnostic classification and evaluation is obtained by identifying specific features 
in order to esctract characteristic patterns which are compared to a library of patterns 
typical of various kinds of heart disorder, and the closeness of the measured signal to 
these patterns. 

5 

Enhanced appreciation and identification of the heart sound features is obtained by 
"placing" the extracted features in a synthetic acoustic environment relying on 
supplying different signals to the ears of a listener by means of headphones. This is 
obtained by means of so-called Head Related Transfer Functions, or HRTF. 

10 

A specific procedure of the invention is characteristic in that first and second heart 
sounds are detected and placed correctly on a time axis by performing the steps of 
feature extraction and classification based on the energy distribution over time in a 
feature time fiinction by the steps of bandpass filtering, followed by instantaneous 
IS power and lowpass filtering of the original phonocardiographic signal. 

An embodiment of the invention is particular in that it comprises the steps of 
extracting the first and second heart sounds by classification according to energy 
levels, eliminating the contribution of the said first and second heart sounds firom the 
20 signal, performing spectral analysis and feature enhancement in order to obtain the 
energy content present in areas of a time-fi-equency representation delimited by 
frequency band times time interval in the form of energy distributions, classifying the 
energy distributions according to pre-defined criteria comparing the energy 
distributions to a catalogue of distributions related to known medical conditions. 

25 

A fiirther advantageousjemhodiment of the. invention for extracting murmur., 
information is particular in that it comprises the steps of obtaining a digital 
representation of heart sound for a predetermined number of seconds, identifying the 
time of occurence of the first and second heart sounds in each cycle, windowing the 
30 parts of heart sounds falling between the first and second heart sounds, and second 
and first heart sounds, respectively; decomposition of the signals into a 
predetermined first number nl of frequency bands, each band being decomposed into 
a predet^muned second number n2 of time-slices, obtaining a systole (SP) and a 
diastole (DP) power vector consisting of the sum of nl powers measured in each of 
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the tCL time slices, for each combination of a frequency band and a time slice, the 
power values from the different systoles being compared, and the median value being 
chosen to be the standard value for a power vector, obtaining a systole (SMF) and a 
diastole (DMF) mean frequency vector by weighting the power value for each of nl 
5 frequency bands with the mean frequency of the corresponding band, summing the 
results and dividing the sum by the corresponding element in the respective systole 
or diastole power vector, while using the time of occurence of the intensity vectors of 
the various classes for classifying the time distribution of murmurs. 

10 A further embodiment of the invention is particular in that it comprises a step 

preceding the step of obtaining systole and diastole murmur intensity vectors SI and 
DI, namely refining the windowing by setting the values of SP, DP, SMF, and DMF 
of the first or last elements equal to the second or last-but-one values, respectively, if 
the values of the first or last elements of the corresponding vectors fulfil 

IS predetermined deviation criteria. 

A further embodiment of the invention is particular in that still further steps are 
included, namely subjecting the signal to double differentiation before 
decomposition, obtaining a systole (SI) and diastole (DI) murmur intensity vector, 
20 respectively, by taking the logarithm of the corresponding SP and DP vectors, 
classifying the obtained logarithmic vectors into murmur intensity classes, and 
comparing the energy distributions to a catalogue of distributions related to known 
medical conditions. 

25 An apparatus for performing the basic procedure of the invention is particular in that 
it confiprises analog -to-dig ital means for converting a heart sound signal into. sampled ... 
data, means for extracting the first and second heart sounds by classification 
according to energy levels, means for eliminating the contribution of the said first 
and second heart sounds from the signal, means for performing spectral analysis, 

30 means for performing feature enhancement, and multiplication means for obtaining 
the energy content present in areas of a time-frequency representation delimited by 
frequency band multiplied by time interval in the form of energy distributions means 
for classifying the energy distributions according to pre-defined criteria, and 
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comparator means for comparing the energy distributions to a catalogue of 
distributions related to known medical conditions. 

An embodiment of the inventive apparatus is particular in that signal processing 
5 means are used to produce a spatial sound distribution beised on frequency, a low 
frequency band being delivered to a first earpiece of a headphone and a high 
frequency band being delivered to a second earpiece of said headphone, the 
frequency bands containing first and second heart sounds and murmur sounds 
respectively. 

10 

A fiirther embodiment of the apparatus is particular in that said signal processing 
means produce a temporal sound distribution, sound signals being first delivered to a 
first earpiece of the headphone and then being delivered to a second earpiece of the 
headphone. 

15 

A fiirther embodiment of the apparatus is particular in that said signal processing 
means comprise at least one Wiener filter. 

^ 

20 The invention will be more fiiUy described in the following with reference to the 
drawing, in which 

Fig. 1 shows a fiinctional block diagram of the complete information extraction 
process according to the invention, 

25 

Fig. 2 shows the steu£^e..pj:QvLding.the first analysis of the heart signal,. 
Fig. 3 shows an original heart signal and its corresponding spectrogram, 
30 Fig. 4 shows the spectrogram of a bandpass filtered heart signal. 

Fig. 5 shows the result of a time marginal distribution of the energy. 



Fig. 6 shows the identification of large and small "hills" in the feature. 
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Fig. 7 shows the identification of SI and S2 irrespective of relative power. 

Figs. 8(a) and 8(b) are representations of the Wiener scenario as applied to synthesis 
5 and analysis, respectively, of filters for right and left ears for stereo headphones 
according to the invention. 

Fig. 9 & 10 are representations, respectively, of a one and a two band approach to the 
application of the Wiener scenario to the stereo headphones according to the 
10 invention. 

Fig. 1 1 is a representation of an alternative manner of presenting heart sounds to 
right and left ears wth the stereo headphones according to the invention, and 

15 Fig. 12 is a modfied embodiment for spatially distributing heart sounds. 

In Fig. 1 is seen a fixnctional block diagram of an embodiment of the procedure and 
sub-procedures according to the invention. The following description relates to a 
practical example of an embodiment according to the invention. 

20 

Example 1. 

The input for the procedure consists of 8 seconds of heart sound signal, sampled at a 
rate of 1000 Hz and read into a digital register subsequent to A/D conversion. The 
procedure is described with reference to modem digital technology, however in 
25 prindple the various classification and sorting of time intervals and levels may be 
performed by analogue means on DC voltages and traditional gates. 

The detector for SI and S2 essentially consists of two separate processes, a feature 
extraction part and a classification part. The purpose of the feature extraction is to 
30 transform the input signal to a domain, in which the respective location in time of SI 
and S2 is more distinct than in the original signal. The classification part determines 
the predse location of SI and 82 and correctly identifies them as such. 



wo 2004/032741 



7 



'CT/DK2003/000680 



10 



15 



20 



25 



In Fig. 2 is demonstrated how murmurs may be observed in the spectrogram of a 
time function of an original heart sound. The spectrogram is obtained by Fast Fourier 
Transform. The first and second heart sounds SI and S2 have only a low-firequency 
content compared to the broad-band nature of the murmurs, and for this reason the 
signal is band-pass filtered by convolution of the original signal with the impulse 
response fimction of a bandpass filter. The corresponding spectrogram is shown in 
Fig. 3, in which peaks of higher energy are visible but not clearly identifiable. In 
order to obtain a time function of the occurence of these higher energies, the time 
marginal distribution of the spectrogram is performed according to Eq. (1): 



Hereby a "final feature" is obtained as a time function as shown in Fig. 5. In essence, 
this time function is obtained by bandpass filtering, instantaneous power extraction 
and lowpass filtering. It is now clear that the "final feature** displays a "hill" every 
time an SI or S2 occurs in the heart signal. 

As the magnitudes of the "hills" corresponding to SI and S2 are comparable, it is 
necessary to distinguish between them by applying classification rules. First all 
"hills" in the "final feature" must be identified. This is obtained for all samples of the 
time function which fulfil the following criteria: 

feature(k-l) < feature (k) and feature(k) > feature (k+1). 

TfteTnexf step is to construct a"faBle of possible systoles. A~sysl61e"is a pair of "hills"* 
(SI and S2) constrained by the time distance between the "hills". The time distance 
must fall within the following limits: 

230 ms < T < 500ms for human hearts. 




The final sequences of systoles is determined by finding the sequence of systoles in 
the table having maximum energy that fulfil the following constraints: 
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- systole time deviation < 18% 

- time between systoles (diastole) > 0.9 times systole time 

- amplitude deviation of Sl< 500% 

- amplitude deviation of S2 < 500% 

5 - in the case of overlapping systoles, the systole with maximum energy must be 
selected. 

The result of the identification is displayed in Fig. 7, in wMch a fat black line to the 
top of a "hill" indicates the time position of a first heart sound SI and a thin black 
10 line a second heart sound S2. 

With the time positions of the first (SI) and second (S2) heart sounds correctty 
detected in the signal (given as sample numbers, corresponding to positions 
measured in milliseconds) it is now possible to evaluate the much weaker sounds, the 
15 heart murmurs. In the following, these detected time positions will be referred to as 
SI markers and S2 nuirkers, respectively. Reference is again made to Fig. 1. 

Delimitation of Systoles and Diastoles. 

Only the systole and diastole parts of the heart sound signal are used for the murmur 
20 detection. All periods, beginning 50 milliseconds after an SI marker and ending 50 
milliseconds before the immediately following S2 marker, are defined as systoles. 

Correspondmgly, all periods, beginning 50 milliseconds after an S2 marker and 
ending 50 milliseconds before tiie immediately following SI marker, are defined as 
25 diastoles. This is a primitive but effiraent manner of eliminating the influence of the 
very . energetic. first.and.secondJieart sounds. At alata-.stagein-the performance of 
the procedure some corrections are made (vide below), but it may be more 
advantageous to perform the elimination using more refined approaches at this early 
stage in tiie procedure. 

30 

Time and Frequency Decomposition of Systoles and Diastoles. 
The sound eneigy content in the sound signal is calculated by means of a 
spectrogram based on the Discrete Fourier Transform (DFT) using a vector length 
which is a power of 2, such as 16. In order to be able to classify murmurs regarding 



wo 2004/032741 




CT/DK2003/000680 



frequency contents and time distribution, each systole and diastole is decomposed 
into 14 frequency bands and 10 time slices, the two lowest frequency bands being 
discarded. The 14 frequency bands cover the frequency range from 62.5 Hz to 500 
Hz, each having a width of 31,25 Hz. 

5 

Befr^re the calculation of the spectrogram, the sound signal is differentiated twice 
(corresponding to a second order high-pass filtration) in order to take into account 
the frequency characteristics of the human hearing, being more sensitive to higher 
than lower frequencies within the frequency range in question. 

10 

It is considered that a parallel bank of band pass filters will perform faster in the 
present environment. 

The 10 time slices for a given systole or diastole all have the same width, 
15 corresponding to 1/10 of the total length of the systole/diastole. 

The combination of frequency bands and time slices creates a 14x10 matrix for each 
systole/diastole. For each element in these matrices, the energy content is divided by 
the width of the relevant time slice, thus yielding matrices containing the heart sound 
20 power (energy per time) for the 140 time/frequency elements of each 
systole/diastole. 

Definition of Standard Systoles and Diastoles for Each Frequency Band. 

The mabices for each systole are combined to a single 14x10 systole (S) matrix by 

25 median filtration: 

For each combination of a frequency range and a time slice, the power values from 
the different systoles are compared, and the median value is chosen to be the 
standard value. This is an efficient way of obtaining a stable value. Thus, for each of 
the 14 frequency bands (rows in the matrix), 10 standard power values combine to a 

30 standard systole. 

The diastole matrices are combined to a D matrix in the same way. 
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Extraction of Power and Frequency Feature Vectors. 

A systole power (SP) vector vnth 10 elements is constructed by summing the 14 
standard power values for each of the 10 time slices. Thus, the SP vector consists of 
the column sums for the S matrix. 

5 

A diastole power vector (DP) is constructed in the same way. 



A systole mean frequency (SMF) vector (also vnth 10 elements) is calculated by 
weighting the power value for each frequency band with the mean frequency of the 
10 corresponding band, summing the 14 results, and dividing the sum with the 
corresponding element in the SP vector. 

Correspondingly, a diastole mean frequency (DMF) vector is calculated. 



15 Correction of Feature Vectors for S1/S2 Remnants. 

Due to the very simple definition of systoles and diastoles, the first and/or last tenths 
of some of the systoles and diastoles may be "contaminated" with parts of SI or S2. 
Typically, this is manifested by larger values of the first/last elements in SP/DP and 
lower values of the corresponding elements in SMF/DMF, because of the high power 
20 and the relatively low frequencies of SI and S2 compared to the murmurs in systoles 
and diastoles. 



Therefore, the beginning and end of the SP, SMF, DP, and DMF vectors are 
examined and corrected if necessary in dependence of the following relationships: 

IfSMF(2)-SMF(l) > 2*|SMF(3)-SMF(2)| or 

SP(1)-SP(2) > 2*1 SP(2) - SP(3) I, 

SP(1) - = SP(2) and SMF(l) = = SMF(2) . 



30 



IfSMF(9)-SMF(10) > 2*1 SMF(8) - SMF(9) I or 

SP(10)-SP(9) > 2*1SP(9)-SP(8)|, 

SP(IO) = = SP(9) and SMF(IO) = = SMF(9) . 
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Corresponding examinations and corrections are performed for DP and DMF. 



Creation of Murmur Intensity Vectors. 

The elements in a systole intensity (SI) vector is created from the elements in the SP 
S vector in the following way using absolute values: 



10 loglO(SP(x)) < -1.25 SI(x) = 0 

-L25 < loglO(SP(x)) < -0.80 SI(x) = 1 

-0.80 < loglO(SP(x)) < -0.35 SI(x) = 2 

15 

-0.35 < logl0(SP(x)) < +0.10 SI(x) = 3 

+0.10 < logl0(SP(x)) < +0.55 SI(x) = 4 

20 +0.55 < loglO(SP(x)) < +1.00 SI(x) = 5 

+1.00 < loglO(SP(x)) SI(x) = 6 

25 A diastole intensity (DI) vector is constructed in the same way. 

It may be relevant to use values relative to e.g. the intensity of SI and/or S2, in 
which case the logarithmic conversion may use other limits than given above. 



30 



Correction of Murmur Intentity Vectors for Noise. 

In order to correct for transient noise signals, the following corrections are 
performed: 
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SI(1) is set to 0, if SI(2) is 0 
SI(IO) is set to 0, if SI(9) is 0 

If any element in SI is more than 1 larger than both of its neighbours, the element is 
S set to be equal to the highest of the neighbours. 

Similar corrections are performed for DI. 

Classification of the Intensity of Murmurs. 

10 The intensities of any systolic and/or diastolic murmur is defined as being the 
maximum value of SI and/or DI, resp. 

If the maximum values are both 0, the heart sound is classified as containing no 
murmurs. 

15 

Classification of the Time Distribution of Murmurs. 

If at least one of the maximum values found is larger than 0, the systolic and/or 
diastolic murmurs are classified according to the profiles of SI and DI, resp. 

20 Any systolic murmur is classified within the first class in the list below whose 
description matches the content of SI: 

Systolic ejection murmurs: The values in SI are increasing to a certain 

point and decreasing after that point. 

25 

Steps that care neither or (^creasing are 
allowed within the increasing as well as the 
decreasing part of the vector. 

30 Early systolic murmurs: The last five values in SI are all 0. 



Early-mid systolic murmurs: 



The last three values in SI are all 0. 
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The first five values in SI are all 0. 



Mid-late systolic murmurs: 



The first three values in SI are all 0. 



5 Pansystolic murmurs: 
descriptions. 



The SI vector does not match any of the above 



Any diastolic murmur is classified within the first class in the list below whose 
10 description matches the content of DI: 

Decreasing diastolic murmurs: The values in DI are decreasing, but never 

increasing. 



15 



Steps that are neither in- or decreasing are 
allowed 



20 



Diastolic murmurs with 
pre-systolic accentuation: 

Uniform diastolic murmurs: 
descriptions. 



DI(9) > DI(8) and DI(10) > DI(9). 



The DI vector does not match any of the above 



25 Classification of the Frequency Range of Murmurs. 

Systolic and diastolic murmur fi'equencies are classified according to the firequency 
band containing the largest power value in the tenth(s) of the systole/diastole 
corresponding to the found maximum values of SI/DI. 



30 



If the largest power value is found in one of the two lowest fi-equency bands 
(containing fi-equencies below 125 Hz), the murmur is classified as a low-fi-equency 
murmur. 
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If the largest power value is found in one of the eight highest frequency bands 
(containing frequencies above 250 Hz), the murmur is classified as a high-frequency 
murmur. 



5 If the none of the above is the case, the murmur is classified as a medium-frequency 
murmur. 



Use of the Murmur Detection output 

The output from the procedure is either a string describing the found murmur(s) or 
10 three values for each found murmur coding for the intensity, the time distribution and 
the frequency range of the munnur(s). The values may either be used for 
classification or form the coordinates in a murmur representation 

It will be noted that in the above description of a specific embodiment that 
15 apparently arbitrary steps were introduced of double differentiation (second-order 
highpass filtering) and of applying a logarithmic frmction in order to obtain intensity 
values. These steps have a psychoacoustic foundation related to the hearing of the 
auscultating person. It is obvious that the classification may well proceed without 
these steps, however it has been determined in practical use that the classifications 
20 obtained by applying these steps are commensurate with observations made by 
trained medical staff, and that the results thereby obtained are much more directly 
applicable to the auscultated phenomena at hand. In this way the medical 
professional will be much fiirther aided than by mere reading and comparing sets of 
three values. 

25 

Example 2. 

The extraction of features may be used in an enhanced manner by creating a 
synthetic spatial environment for two ears to listen to via headphones. The ability to 
distinguish more clearly between several simultaneously occuring phenomena has 
30 been used by jet airplane pilots to increase separation when listening to several 
communication channels simultaneously. The use of special filters for these 
purposes is known and research has been conducted to determine the filter effect of 
the head (characterized as Head related Transfer Functions, HRTFs) with regard to 
sound from two microphones based on the distance and source; see, HRTF 
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15 



Measurements of a KEMAR Dummy-Head Microphone, Bill Gardner and Keith 
Martin, MIT Media Lab, (HTTP://sound.media.mit.edu/BCEMAR/hrtfdoc.txt, August 
8, 2000). However, to date, no known attempt has been made to use this ability in 
connection with features extracted from a heart signal and presented as index 
5 markers simultaneously with the heart signal itself so as to present different features 
to each of the physician's ears. 

Figures 8(a) & 8(b) show one manner in which filters can be used to spatially 
distribute the sounds delivered to the physician's ears to reflect the angular difference 
10 between sounds received by the left and right ears. In particular, either the signal to 
both ears can be filtered or the original signal can be preserved and sent to one ear. 
To do so, it is necessary to synthesize the difference at each angle, and this can be 
done by applying the known Wiener scenario, by which an estimate of the optimal 
filter coeflBcients for a specified FIR filter can be arrived at. 

15 

In Figures 8(a) & 8(b), destimaie(n) designates the impulse response signal z"^ to be 
synthesized for the ear farthest away from the source, and u(n) designates the 
impulse response for the other ear, which by proper filtering should mimic destimatc(n). 
By means of a least mean square algorithm, the coefficients w are adjusted to 

20 minimize the error signal e(n). These special filters, one for each angle, are used to 
arbitrarily place a sound spatially when listening via a stereo headset, the original 
sound being presented to one ear and a filtered version to the other ear. 
The listening experience can also be expanded by creation of a synthetic listening 
space in which low frequency sounds, such as heart beats with the enhancements 

25 described in Example 1, are perceived as coming from, e.g., the left side while high 
frequency sounds, such as heart murmurs, from- the right side, -Similarly, eadier 
phenomena could be made to appear on, e.g., the left side and subsequent phenomena 
on the right so that, with a repetitive sequence, there would be a repetition of sounds 
moving from left to right. In these manners, separating and distinguishing of features 

30 is facilitated. Figures 9 & 10 represent one band and two band scenarios, 

respectively, for achieving these effects. In Figure 9, the input sound passes through 
to the left ear, while the signal to the right ear is processed in one of the above 
manners. In Figure 10, the input sounds are separated and independently processed 
en route to each ear. 
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Figure 11 shows an arrangement for transformation of heart sounds from frequency 
distribution to a spatial distribution. The sound signal is first divided into a number 
of frequency bands by normal filters or orthogonal filters, orthogonal filters 

5 preventing redundancy, which ensures energy preservation. The output from each 
filter has a direct path and a delayed path to the matrix circuits for the left and right 
channels. In the matrices, a weighted sum of the input signals is formed in such a 
way that the lowest to highest frequency bands are perceived as being spatially 
distributed from left to right, when played back via a stereo speaker system or 

10 preferably via a stereo headset. In this way, an alternative presentation can be 
offered which adds a new dimension to the sound, which apparently enhances the 
perceived frequency resolution, and by that, the ability to recognize murmurs etc. 

Figure 12, on the other hand, shows an arrangement for the transformation heart 
15 sounds from a temporal distribution to a spatial distribution. In this embodiment, the 
spatial location of the sound follows the temporal location in the heart cycle, from a 
first heart sound to the next first heart sound, si to si. This means that the systole, si 
to s2, is located on the left side and the diastole, s2 to si, on the right side. It is like 
the balance is automatically adjusted with time triggered by a signal derived from the 
20 input signal itself. In this way, physicians are offered an alternative presentation 
which is meant to help locate murmurs in the heart cycle, mainly systolic or mainly 
diastolic. 



It will be understood that once the signal has been converted to digital representation 
25 of data, its manipulation may take place in dedicated processors, RISC processors, or 
general purpose computers,-the outcome-of thejtianipulation being-solely dependent 
on the instructions performed on the data under the control of the program written for 
the processor in order to obtain the fimction. The physical location of the data at any 
one instant (i.e. in varying degrees of processing) may or not be related to a 
30 particular block in the block diagram, but the representation of the invention in the 
form of interconnected fijnctional blocks provides the skilled person with sufficient 
information to obtain the advantages of the invention. 



wo 2004/032741 ^VCT/DK2003/000680 

17 

The foregoing description of the specific embodiments will so folly reveal the 
general nature of the present invention that others skilled in the art can, by applying 
current knowledge, readily modify or adapt for various applications such specific 
embodiments without undue experimentation and without departing fi-om the generic 

5 concept, and therefore, such adaptations and modifications should and are intended 
to be comprehended within the meaning and range of equivalents of the disclosed 
embodiments. It is to be understood that the phraseology or terminology employed 
herein is for the purpose of description and not of limitation. The means, materials, 
and steps for canying out various disclosed fonctions may take a variety of forms 

10 without departing firom the invention. 



Thus, the expressions "means to ... " and "means for ...", or any method step 
language, as may be found in the specification above and/or in the cldms below, 
followed by a fonctional statement, are intended to define and cover whatever 

15 structural, physical, chemical, or electrical element or structure, or whatever method 
step, which may now or in the future exist which carries out the recited fonctions, 
whether or not precisely equivalent to the embodiment or embodiments disclosed in 
the ^ecification above, i.e., other means or steps for carrying out the same fonction 
can be used; and it is intended that such expressions be given their broadest 

20 interpretation. 
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PATENT CLAIMS 

1. A procedure for extracting information from a phonocardiographic signal 
obtained from a transducer and subjected to signal processing, 

5 characterised in that first and second heart sounds are detected and placed 
correctly on a time axis by performing the steps of feature extraction and 
classification based on the energy distribution over time in a feature time function by 
the steps of bandpass filtering, followed by instantaneous power and lowpass 
filtering of the original phonocardiographic signal. 

10 

2. A procedure for extracting information from a phonocardiographic signal 
according to claim 1, including identification of characteristic signal components, 
characterised in that it comprises the following steps: 

- extracting the first and second heart sounds by classification according to energy 
15 levels, 

- eliminating the contribution of the said first and second heart sounds from 
the signal, 

- performing spectral analysis and feature enhancement in order to obtain the energy 
content present in areas of a time-frequency representation delimited by frequency 

20 band times time interval in the form of energy distributions 

- classifying the energy distributions according to pre-defined criteria 

- comparing the energy distributions to a catalogue of distributions related to known 
medical conditions. 



25 3. A procedure used for extracting murmur information, 

..character! s.e.d_.In that it comprises the-fi)llowing steps: 

- obtaining a digital representation of heart sound for a predetermined number of 
seconds, 

- identifying the time of occurence of the first and second heart sounds in each cycle, 
30 - windowing the parts of heart sounds falling between the first and second heart 

sounds, and second and first heart sounds, respectively 

- decomposition of the signals into a predetermined first number nl of frequency 
bands, each band being decomposed into a predetermined second number n2 of 
time-slices 
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- obtaining a systole (SP) and a diastole (DP) power vector consisting of the sum of 
nl powers measured in each of the n2 time slices 

- for each combination of a frequency band and a time slice, the power values from 
the different systoles are compared, and the median value is chosen to be the 

5 standard value for a power vector 

- obtaining a systole (SMF) and a diastole (DMF) mean frequency vector by 
weighting the power value for each of nl frequency bands with the mean frequency 
of the corresponding band, summing the results and dividing the sum by the 
corresponding element in the respective systole or diastole power vector 

10 - while using the time of occurence of the intensity vectors of the various classes for 
classifying the time distribution of murmurs. 

4. A procedure for extracting murmur information according to claim 3, 
characterised in that a step preceding the step of obtsuning systole and 
15 diastole murmur intensity vectors SI and DI consists of refining the windowing by 
setting the values of SP, DP, SMF, and DMF of the first or last elements equal to the 
second or last-but-one values, respectively, if the values of the first or last elements 
of the corresponding vectors fulfil predetermined deviation criteria. 



20 5. A procedure according to claim 3, characterised in thatfiirther 

steps are included in the procedure, comprising 

- subjecting the signal to double differentiation before decomposition 

- obtaining a systole (SI) and diastole (DI) murmur intensity vector, respectively, by 
taking the logarithm of the corresponding SP and DP vectors, 

25 - classifying the obtained logarithmic vectors into murmur intensity classes 

-jcomparing the energy-distributionsJo a catalogue of distributionsjrelatedjtaknown. 
medical conditions. 



6. An apparatus for performing the procedure according to claim 1, 
30 characterised in that it comprises analog-to-digital means for converting a 
heart sound signal into sampled data, means for extracting the first and second heart 
sounds by classification according to energy levels, means for eliminating the 
contribution of the said first and second heart sounds from the signal, means for 
performing spectral analysis, means for performing feature enhancement, and 
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multiplication means for obtaining the energy content present in areas of a time- 
frequency representation delimited by frequency band multiplied by time interval in 
the form of energy distributions means for classifying the energy distributions 
according to pre-defined criteria, and comparator means for comparing the energy 
5 distributions to a catalogue of distributions related to known medical conditions. 



7. An apparatus for performing the procedure according to claim 1, wherein 
signal processing means are used to produce a spatial sound distribution based on 
frequency, a low frequency band being delivered to a first earpiece of a headphone 
10 and a high frequency band being delivered to a second earpiece of said headphone, 
the frequency bands containing first and second heart sounds and murmur sounds 
respectively. 



8. An apparatus according to claim 6, wherein said signal processing means 
15 produce a temporal sound distribution, sound signals being first delivered to a first 
earpiece of the headphone and then being delivered to a second earpiece of the 
headphone. 



20 



9. A apparatus according to claim 7 or 8, wherein said signal processing 
means comprise at least one Wiener filter. 
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